-3 - 
Amendments 



Frank GUARNIERI 
Appl. No. 09/183,267 



In the Specification: 

On page 4, on the line following the heading "Brief Description of the Drawings," 
please add the following paragraph: 

The patent or apphcation file contains at least one drawing executed in color. 
Copies of this patent or patent application publication with color drawing(s) will be 
provided by the Office upon request and payment of the necessary fee. 

Please replace the paragraph beginning on page 12, line 20 with the following paragraph: 
hi Figures 5 A and 5B, solutions obtained with co-crystals of elastase inhibitors are 
compared with data obtained by the methods herein described. Li Figure 5 A, the solutions 
for six co-crystallized inhibitors are shown, with the inhibitor molecules overlaid on each 
other (non-space-filling representation, with the elastase segment represented by a space- 
filling illustration). These inhibitors are trifluoroacetyl-l-lysyl-l-prolyl-p-isopropylanilide 
( Mattos, C. etal.Nat Struct Biol 7:55-58 ri994): crystal solution: PDBID:1ELA Mattos 
a al., as submi t ted Apiil 30, 1994 ), trifluoroacetyl-l-lysyl-l-leucyl-p-isopropylanilide 
f Mattos, C. etaL.Nat. Struct Biol 7:55-58 (1994): crystal solution: PDB ID: lELB Mattos 
a al., as submitted June 22, 1 9 94 ), trifluoroacetyl-l-phenylalanyl-p-isopropylanilide 
( Mattos. C. etaL.Nat Struct Biol 7:55-58 (1994): crystal solution: PDB ID: lELC Mattos 
ct — ah; — as — submitted — April — 30; — 1:994), trifluoroacetyl-l-phenylalanyl-l-alanyl- 
p-trif1iinrnmethy1aninide p-trifluoromethvlanilide ( MdXtos.C. et al. Biochemistr y 3 4:319 3- 
203(1995): crystal solution: PDB ID: lELD Mattos ctaL, as submitted Fcbiuai-y 14, 1 99 5 ), 



- 4 - Frank GUARNIERI 

Appl No. 09/183,267 

trifluoroacetyl-l-valyl-l-alanyl-p-trifluoromethylanilide ( Mattos, C. et aL. 
Biochemistry i-/:3193-203 (1995^: crystal solution: PDB ID: lELE Mattos ct ah, as 
submitted February 14, 1995 ) and n-(tert-butoxycarbonyl-alanyl-alanyl)-o-(p-nitrobenzoyl) 
hydroxylamine (Ding, X. et al. Biochemistry 3^:7749-56 ( 1 995): crystal solution: PDB ID: 
lELF Ding ct al., as submitted July 10, 1995 ). In Figure 5B, the solutions for approximately 
10 ORFs, which are in their respective high affinity protein binding states are overlaid. 
Both methods identify a region which favors the binding of aromatic moieties. The 
simulation process achieves approximately 90% 3D geometric identity with the 
crystallography results. 

Please replace the paragraph beginning on page 17, line 2 with the following paragraph: 
Aspects of the simulations used in the invention can be illustrated with calculations 
used to determine the strength of water binding to a synthetic polynucleotide (Guamieri, F.. 
and Mezei. M.. J. Am. Chem. Soc, 77^:8493-8494 (1996)) .rrMl This illustration can be 
described as follows: 

Please replace the paragraph beginning on page 17, line 10 with the following paragraph: 
Grand canonical ensemble simulations are generally performed by placing a 
molecule in a periodic simulation cell, setting a parameter B, which is representative of fi-ee 
energy, in such a way as to achieve an experimentally determined density, sampling 
potential hydration positions around the molecule by inserting and deleting water molecules 
from the simulation cell using a technique such as cavity-bias (Mezei. M.. Mol. Phys. 
57:565-582 (1994^: Resat. H.. and Mezei, M., J. Am. Chem. Soc. 776:7451-7452 
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C1994}},[[^'^]] and accepting or rejecting the attempt based on a Metropolis Monte Carlo 
(Metropolis, N. et aL. J, Chem, Phvs. 27:1087-1092 ri9531X m criteria using a grand 
canonical ensemble probability function (Tolman. R„ in The Principles of Statistical 
Mechanics. Dover Press, New York (1971)) . [[^]] The parameter B is related to the excess 
chemical potential ju' as follows: B = /u'/kT + ln<N>, where k is Boltzmann*s constant, T 
is the absolute temperature, and <A^ is the mean number of molecules of the ORF, which 
here is H2O. In the method of simulated annealing of chemical potential, the simulation is 
started with a large initial 5- value so that a higher percentage of water insertion attempts 
are accepted. This causes the simulation cell to be flooded with water molecules. After this 
grand canonical ensemble simulation at high excess chemical potential is equilibrated, 
subsequent simulations are carried out at successively lower 5-values. This successive 
lowering of the 5-values causes a gradual removal of the bulk water molecules from the 
simulation cell. As the chemical potential is fiirther "annealed", a point is reached at which 
water molecules do not readily leave the cell, thereby identifying those water molecules that 
are strongly influenced by the DNA, the so-called "bound water molecules". As the excess 
chemical potential is again lowered, ultimately some of these bound waters start to leave the 
cell. Since chemical potential is a free energy, this simulated annealing of chemical 
potential yields a numerical estimate of the differential free energy of binding of the 
different bound water molecules. It must be emphasized that our utiUzation of the term 
"annealing" applies strictly to the value of the chemical potential and that the temperature 
is kept constant at, for example, 298 K in all the simulations. For all simulations the DNA 
was held fixed, water molecules were added and deleted throughout all parts of the cell, 
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extensive canonical Monte Carlo was performed between accepted grand canonical Monte 
Carlo steps, and periodic boundary conditions were used. 

Please replace the paragraph beginning on page 18, line 6 with the following paragraph: 

As an illustration of the method, a simulated annealing of chemical potential on a 
d(CGCGAATTCGCG)2 was performed, starting with 5 = 1.0 down to -26 in 37 increments 
performing 2,000,000 cavity-biased grand canonical ensemble Monte Carlo steps at each 
5-value. The final configuration of the simulation with B = -6, has 1120 boxmd water 
molecules. The final configuration of the simulation with B = -8, has 533 bound water 
molecules. The final configuration of the simulation with B = -9, has 390 bound water 
molecules. The final configuration of the simulation with 5 = -1 1, has 215 bound water 
molecules. The most salient feature of this progression is the differential hydration of the 
major and minor groove of the DNA. The B = -6 simulation shows the DNA essentially 
uniformly solvated. The J? = -8 simulation clearly shows that upon lowering of the chemical 
potential by 2 B-units, a majority of the nonbulk extracted waters come from the major 
groove, while the minor groove remains almost unaffected. Annealing the chemical 
potential fiirther (B = -9) still leaves the minor groove well hydrated while the major groove 
is almost stripped. Lowering B even fiirther (B = Al) results in the removal of almost all 
water molecules from both the major and minor groove. Quantitation of the hydration of 
the DNA as a fimction of chemical potential was computed by proximity analysis 
rMehrotra.P.K.,andBeveridee.D.L., J.^m. Chem. Soc, 702:4287 (1980^: Mezei, M.,Mo/. 
SimuL 7:327-332 (1988) (describing the effects of different partial charges on proximitv 
analysis)) !" ["^'^11 with the results shown in Table 1: 
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first hydration shell first and second hydration shell 

minor groove ma.ior groove minor groove maior groove 

no. of no. of no. of no. of 



B 


waters 


density 


waters 


density 


waters 


density 


waters 


density 


-6 


7.27 


0.013 


13.23 


0.012 


21.3 


0.021 


41.7 


0.011 


-8 


5.4 


0.010 


5.06 


0.004 


14.6 


0.015 


11.8 


0.003 


-9 


4.08 


0.007 


4.36 


0.004 


11.5 


0.011 


9.7 


0.003 


-11 


1.04 


0.002 


2.11 


0.002 


3.9 


0.004 


4.2 


0.001 



For B = -6, the first hydration shell (defined by the position of the first minimum of the 
radial distribution function) of the major and minor groove has a comparable density (0.01 2 
and 0.013, respectively), while the second hydration shell of the minor groove has twice the 
density of the major groove. ForB = -8 the hydration difference becomes quite pronovmced 
with the minor groove first and second shell hydration density being 2.5 fold and 5 fold 
higher than the major groove, respectively. For 5 = -11 the major and minor groove 
hydration density again becomes equal because at this value of the excess chemical potential 
both grooves are essentially stripped bare. 

Please replace the paragraph beginning on page 19, line 6 with the following paragraph: 

Illustrating the differential hydration propensities of the major and minor grooves 
of DNA is computationally imdemanding (3 days of CPU time to run one annealing 
schedule and 3 days of CPU time to run one proximity analysis (Mezei, M.. and Beveridge, 
DX., in Methods in Enzvmolosv. Packer, ed.. Academic Press, New York, pp. 21-47 (1986) 
(describing the effects of volume element calculations — which can be CPU intensive — on 
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proximitY analvsis')) [[^]] on an SGI Power Challenge) using simulated annealing of 
chemical potential because only a coarse "cooling" schedule of the chemical potential is 
required. Since the chemical potential is a free energy, a very fine cooling schedule may be 
used to estimate quantitatively the hydration free energy difference of two different 
functional groups or even two different atoms of the DNA. Two atoms that desolvate at the 
same 5-value have similar solvation free energy, or altematively, require a finer cooling 
schedule to resolve the differences. It should be noted that the model system used here 
consisted of ionic DNA with 22 negative charges and no sodium counterions. The findings 
presented herein about the preferential hydration of the minor groove corresponds very well 
to results from X-ray crystallographic and NMR studies. Possible reasons for the stronger 
binding of water molecules in the minor groove may include the following: the high density 
of the charged rows of phosphate groups, steric constraints, and specific water — water, 
water — DNA interactions. 

Please replace the paragraph beginning on page 21, line 6 with the following paragraph: 
References 

All publications and references, including but not limited to patents and patent 
applications, cited in this specification are herein incorporated by reference in their entirety 
as if each individual publication or reference were specifically and individually indicated 
to be incorporated by reference herein as being fiiUy set forth. Any patent application to 
which this application claims priority is also incorporated by reference herein in its entirety 
in the manner described above for publications and references. The cited documents 
incorporated by reference into this disclosur e include : 
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